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Abstract 

We discuss the possibility of deciding on the symmetry of the superconducting 
phase in the organic superconductors (Bechgaard salts), using tunnelling experiments. 
We first briefly review the properties of organic superconductors, and the possibility 
to have unconventional (triplet) superconductivity in these systems. We then present 
a simple scheme for computing the full current-voltage characteristics for tunnelling 
experiments within the framework of the non-equilibrium Keldysh Green function for- 
malism. This formalism is flexible enough to address different pairing symmetries 
combined with magnetic fields and finite temperatures at arbitrary bias voltages. We 
then discuss extensively how to apply these results to probe for the symmetry of the 
superconducting order parameter. 
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§1. Introduction 



Organic materials have provided physicists with an extraordinary laboratory to study 
the effects of interactions in solids. These materials indeed show remarkable effects due to 
the interplay of interactions and dimensionality. ^ Indeed such compounds made of cou- 
pled chains allow to realize practically the physics of Luttinger liquids, 2 ) one of the very 
few controlled examples of non-Fermi liquids 3 ) due to interactions. But contrarily to other 
realizations of one dimensional systems such as nanotubes, 4 ) -6 ) quantum wires 7 ) -9 ) or edge 
states in the quantum Hall effect, 10 H 2 ) the organics offer unique challenges. Indeed, because 
of their very three dimensional nature, they provide not a single one dimensional electron 
gas, but a very large number of such one dimensional systems coupled together. This allows 
thus for a unique new physics to emerge where the system is able to crossover from a one 
dimensional behavior to a more conventional three dimensional one. 

As a consequence of this dimensional crossover, at low temperature these materials un- 
dergo instabilities towards three-dimensionally ordered states, such as spin-Peierls, antiferro- 
magnetic and even superconducting states. Needless to say, the presence of superconductivity 
in these compounds is a tantalizing and challenging question. Despite a period of quarter of 
a century since the discovery of superconductivity in these materials, 13 ) the mechanism and 
even the symmetry of this superconducting phase have remained elusive, and many efforts 
have been devoted to this subject. Recently, the case for triplet nature of this supercon- 
ducting phase has been made, 14 ) mostly by measurements of the upper critical field and by 
NMR measurements, but the subject is far from being settled. If true, this would put the 
organics in the relatively small club of condensed matter materials where interactions can 
lead to unusual triplet superconductivity, 15 ) other examples being 3 He (a superfluid) and 
Sr 2 Ru0 4 , 16 )' 17 ) and can certainly shed light on the paring mechanism. 

In these notes we explain how point contact tunnelling experiments can be used to probe 
for the symmetry of a superconducting state. In the recent years, the possibilities to perform 
point contact tunnelling have been drastically enhanced with the development of scanning 
tunnelling microscopy (STM). 18 )' 19 ) Correspondingly, theories to interpret tunnelling exper- 
iments in superconductors have evolved from simple semiconducting band models, 20 )' 21 ) to 
more systematic approaches based on the tunnelling Hamiltonian. 22 ) -24 ) We show here how 
one can extend and simplify the formalism to make it more versatile and easy to imple- 
ment. 25 )' 26 ) This allows to study the case of superconducting (singlet or triplet) leads, as 
well as the effects of magnetic fields on the junction. The resulting theory can be thus di- 
rectly used as a probe of the symmetry of the leads. We thus discuss how such a probe can 
be used for the case of the organics. 
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The plan of these notes is as follows. In Sec. El we briefly recall the salient points of 
the physics of organic compounds. In Sec. H3 we give an outline of our formalism for the 
tunnelling between triplet superconductors. 25 )> 26 ) Sec. |U discusses how to test for triplet 
superconductivity in the organics. Finally some conclusions and perspectives can be found 
in Sec. El 

§2. A brief word on organics 

Let us briefly recall some of the properties of the organic materials. We will focuss 
here on the Bechgaard salts TMTSF 2 X which were the first organic compounds to exhibit 
superconductivity, and have thus been the focus of intense experimental and theoretical 
studies. A general review on these, and parent compounds (so called Fabre salts TMTTF 2 X), 
can be found in Ref. 1). 

In addition to the superconducting phase itself, these materials have a notoriously rich 
phase diagram (cf. Fig. and exhibit a host of remarkable properties (non-FL metallic 
behavior, quantized Hall conductance, Frohlich conductivity), many of which are still poorly 
understood. These materials exhibit a quasi-one dimensional structure due to the nature 
of the molecular orbitals protruding from the flat molecule TMTSF (tetramethyltetrasele- 
nafulvalene) which is the basic building block of the Bechgaard salts: the overlap of the 
7i- orbitals of the selenium or sulfur atoms leads to a high mobility of electrons along the 
stacking direction. The hopping integrals in the perpendicular directions are smaller by 
more than one order of magnitude. Estimated values of the hopping integrals along the 
stack direction (a-axis) and the two perpendicular axes pointing towards neighboring stacks 
(b-axis) and towards the anions (c-axis) respectively are: t a : t& : t c — lOOOfT : 100i^ : 30K. 
Therefore one can think of these materials as one-dimensional chains coupled by small inter- 
chain hoppings. Given the hierarchy of transverse coupling the system is first expected to 
become two dimensional and then three dimensional at low temperatures. The normal phase 
of these materials presents unusual transport 28 ^ 30 -* and optical conductivity properties 31 **' 32 ) 
that have been interpreted as the signature of Luttinger liquid physics. 32 )' 33 ) The high tem- 
perature phase of these materials thus clearly shows effects of strong correlations. As the 
temperature is lowered, the material recovers features that are more and more reminiscent of 
a normal Fermi-liquid material, and finally undergoes instabilities towards various ordered 
phases (spin-Peierls (SP), antiferromagnetic (AF), spin-density wave(SDW)) and supercon- 
ducting (SC)). The nature of the molecule (TMTTF vs TMTSF) or the ions, as well as the 
application of pressure allows to control the phase diagram, by changing the hopping inte- 
grals. This modifies the relative importance of the kinetic energy and Coulomb interaction. 
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Fig. 1. Unified experimental phase diagram for the TM compounds (from Ref. 27)). Either pres- 
sure or chemical changes (increasing pressure corresponds to either going from the TMTTF to 
the TMTSF family or changing the anions) yields the same phases [MI: Mott insulator, LL: 
Luttinger liquid metal, FL: Fermi liquid metal, SP: spin-Peierls, AF: antiferromagnetic spin- 
density wave, SC: superconducting]. Note that the pressure axis does not change the doping of 
the material that remains a quarter filled system. The TMTTF family is insulating at ambi- 
ent pressure whereas the TMTSF family shows good metallic behavior at room temperature. 
The superconducting phase, stemming out of the antiferromagnetic one is still very poorly 
understood. 



The chemical and pressure changes have similar effects, which can be summarized by the 
unified phase diagram of Fig. ^ More issues on the effect of interactions and the properties 
of the normal phases can be found in Ref. 34) and references therein. 

Among the ordered phases, the superconducting one 13 -* is by far the most mysterious. The 
mechanism behind it is still unknown, and considerable debate takes place on that point (see 
e.g. Ref. 35)-37) and references therein). Since at the temperature when superconductivity 
occurs, the compound is not in the one dimensional limit anymore, but well in the three 
dimensional regime, one could in principle think of a conventional BCS mechanism. However 
the proximity of the AF phase as well as the fact, clear from the high temperature phase, 
that interactions are particularly strong in these compounds make such a simple pairing 
mechanism very unlikely. Besides the mechanism itself, the symmetry of its superconducting 
order parameter seems to be spin-triplet ^-wave 38 -* but claims of /-wave orbital symmetry 
were also put forward; 39 ^' 40 ) the microscopic origin for the pairing remains largely a matter of 
debate. One of the most striking features is the existence of a common boundary between an 
antiferromagnetic phase (SDW) and the superconducting phase that is stabilized as pressure 
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is increased. Remarkably, the largest values of T c are actually reached near this boundary 
at which the system is expected to display an enhanced SO (4) symmetry. 41 ) Intriguingly, 
recent experiments have identified a pressure window around the boundary at which the 
SDW and the superconducting regions seem to segregate. 42 ) Among the properties of the 
superconducting phase that would be in favor of triplet superconductivity are the sensitivity 
to non- magnetic impurities,, 37 -* the absence of Knight shift signals at the superconducting 
transitions 43 -* and, perhaps the most striking, the anomalously large values of the measured 
upper critical fields. 44 -* Given the nature of these systems, it is difficult to perform phase 
sensitive experiments such as the ones done for the cuprates, so other probes are needed to 
unambiguously decide on the symmetry of the order parameter. We now examine tunnelling 
as such a probe. 



The earliest successful theoretical models to study superconducting tunnelling junctions 
were based on a scattering picture and semiconducting-like bands. 20 )' 21 )' 45 )' 46 ) It was later 
shown that those results can be recovered using a tunnelling Hamiltonian as the starting 
point. 22 )~ 24 ) A large series of experiments on atomic-size contacts 47 ) -52 ) showed impressive 
agreement with the theory by considering a small number of independent conduction chan- 
nels, each of them well described by a point-contact model. Tunnelling is thus a very efficient 
probe of the properties of the junctions, and one can expect to use it to determine the sym- 
metry of the superconducting order parameter in the junction leads. However the techniques 
used so far are either semi-phenomenological or very heavy, and one thus need to find a for- 
mulation that is both simple enough to be extended to the case of unusual superconductors 
and to finite magnetic fields and finite temperature, and yet accurate enough to compute 
the full current-voltage curve. We describe here such a technique, based on the Keldysh 
formalism. More details can be found in Ref. 25), 26). 

Let us use a tunnelling Hamiltonian formalism to calculate the full current- volt age char- 
acteristics of different types of tunnel junctions where each side of the junction can be either 
a normal metal (N), a singlet (S) or a triplet (T) superconductor. We model the system 
with a (one-dimensional) Hamiltonian that includes the two leads and a tunnelling term: 
H = Hi + H 2 + H tun . Each lead is described by 



where K = H — fiN and /i is the corresponding electrochemical potential. All the indices 
are summed over, in particular k is the lattice momentum, Greek indices correspond to the 



§3. Tunnelling 




(3-1) 
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spin, and c G (L, R) = (— 1, +1) sums over the two possible chiralities. £ c fc CT = cvpk — fi — ah 
are the corresponding linear dispersions, shifted by the inclusion of chemical potential and 
magnetic field along the z-axis (for convenience we will take vp = 1). This is the extension 
of the pairing-approximation Hamiltonian found in BCS theory to the triplet case (for more 
details see Refs. 25), 26)). The third term in the Hamiltonian describes the tunnelling: 

#tun = *l (°) (°) • ( 3 ' 2 ) 

i,t',C 

According to this term an electron with spin a can hop from lead £' into lead £ with a 
tunnelling matrix amplitud given by 

t K - ^ Q j . (3-3) 

Since the number of particles in each lead is a conserved quantity in the absence of tunnelling 
(pair fluctuations are to be regarded as an artifact of the mean field approach), we can define 
the current as given by the rate of change in the relative particle number caused by tunnelling. 
One writes 

J = ^([/7 tun ,iV 1 -iV 2 ]) . (34) 

Having in mind that very simple models of the leads suffice to achieve even quantitative 
agreement with the experiment when calculations are carried out to capture the main features 
of point-contact transport on conventional superconductors, with dimensionality playing 
little or no role, it is justified to use one- dimensional leads to carry out all the standard 
calculations. For unconventional superconductors the situation is more complex, because the 
anisotropic nature of the pair wave-function has to be taken into account when modelling the 
leads. But the organic superconductors that we are interested in are supposed to have p-wave 
symmetry and be highly anisotropic. We can therefore, as a first approximation, conveniently 
restrict ourselves to a one-dimensional model and adopt a formalism that encompasses both 
s- and p-wave symmetries, as well as the normal state. 

In order to deal with an out of equilibrium situation, we use the so-called Keldysh for- 
malism. 53 ) We treat the tunnelling term to all orders to calculate the full I-V line and give 
a quantitative account of the subgap structure. Notice that in this framework the current is 
given by 

1 = | E / % (0) VV (0) - < (0) VV (0)) kel (3-5) 

where 'kel' denotes the Keldysh component of the correlation function. Since the current 
depends only on the fields at x = one can integrate the x dependence in the leads to 
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obtain from H a local quadratic action for the contacts: Sj un = Si + S2 + S tU n- With S tU n 
obtained directly from H tnn and Si = j |^\E , J(u;)g~ 1 \E , e (uj). Here is an 8 component 
spinor and g~ l is a matrix whose inverse (g) is given by the standard advanced, retarded 
and Keldysh components of the local Green functions of the lead. For instance, for the case 
where A\ = A 2 = 0, 



[ret,adv] 

"c<T,C(7 



[ret,adv] 



— (uj — fi£ + cah ± i0~ 



2\/\A + caA 3 \ z - (u- M + cah±i0+y 
(A + caA 3 t ] ^ 



2\ \A + caA 3 \ - (u - ^ + cah ± iO+Y 



And the Keldysh component is g kcl = (g TCt — g adv ) tanh ((cu — fii) /2T). 

Attention must be paid to the fact that frequencies have different reference levels in each 
side of the junction when a bias is present. While, in each side of the junction, states of 
frequencies with equal positive and negative shifts from the Fermi level are related by the 
pairing fluctuations, across the junction, same frequency states are related by tunnelling (see 
Fig. |2J). Thus, the full action for the junction is not diagonal. To each value (u ) in the 




Fig. 2. Representation of a set of frequencies involved in a high-order tunnelling process (the 
vertical axis corresponds to frequencies). The chemical potentials, thermal distributions and 
superconducting quasiparticle densities of states in each side of the junction are all indicated. 
The figure illustrates the lowest order contribution to tunnelling at low voltages for a case that 
corresponds to a process involving one electron plus three Cooper pairs. 



frequency window (of size eV — H\ — ^2) defined by the chemical potentials in the two leads, 
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one infinite set of related frequencies can be assigned (p > 0): 

LOp = 2/i2-pmod2 — 

The action is block diagonal between these sets. Discretizing the frequencies in the window 
defined by the voltage difference defines in turn a discretization of the whole frequency axis. 
We shall deal with one set of frequencies at a time, and since the sets are infinite, we will 
truncate their hierarchy. This amounts to introducing a limit in the number of Andreev 
reflections. One can then numerically invert the action block by block (for more details 
see Refs. 25), 26)). The off-diagonal Green function matrix elements thus obtained allow to 
compute the current using Eq. ()3-5|) . 

§4. Test for triplet 

As it was stated above, in the Bechgaard salts, both the orbital and the spin symmetries of 
the superconducting phase are not clearly known, not to speak of the microscopic mechanism 
responsible for the pairing. Since the uncertainties are such, we think the best way to proceed 
would be if the experiments try to resolve these issues one at a time and to our opinion the 
most approachable one, from the point of view of tunnelling, is the issue of the symmetry 
in spin space of the superconducting order parameter. Thus, the question to be answered is 
whether the electrons in a Cooper pair form a singlet or a triplet. 

To discuss and compare the I-V characteristics for different types of junctions, we choose 
some convenient set of parameters that we will use in all the figures. For the tunnelling 
overlap integral we choose the values t = 0.2 and t = 0.5 (that correspond, in the notation 
of Ref. 20), to a ~ 0.15 or Z = 2.4 and to a = 0.64 or Z = 0.75, respectively), and 
when there is a magnetic field we fix its value to h = 0.2 in units of A (by A we mean 
the magnitude of the singlet gap, Z\ , or of the triplet vector order parameter depending on 
the case). We show curves for the dc response in the limit of vanishing temperatures. The 
geometry of the junction we consider corresponds to tunnelling perpendicular to the chains 
of the quasi-onedimensional compound, with an orbital order parameter aligned along them. 
In this situation no mid-gap states are expected and, therefore, no midgap features either. 
We concentrate on what happens at the conduction edge, in particular the effects of applied 
fields. We refer the reader to the literature for some recent studies that look at the effect of 
fields on zero bias anomalies. 54 ) 

Let us start with the case of normal- metal-superconductor junctions, that would core- 
spond to standard STM experiments. We show in Fig.EJ-(a) typical curves for a point-contact 
junction between a normal metal and a conventional singlet-paring superconductor. The 
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Fig. 3. I-V characteristics of normal-superconductor junctions, (a) N-S junctions for t = 0.2 (lower 
curves: with and without applied magnetic field, dashed and solid line respectively, h = 0.2) 
and t = 0.5 (upper curve, solid line only); the curves are vertically displaced for clarity, (b) 
N-T junctions for t = 0.2 (upper curves) and t = 0.5 (lower curve). 

solid lines correspond to the N-S junction in zero field and the dashed line is for one of the 
junctions in the presence of a magnetic field that produces what would be seen as a Zeeman 
splitting of the differential conductance peak. The second part of the figure corresponds to a 
junction between a normal metal and an unconventional triplet-pairing superconductor. The 
solid lines correspond to the N-T junction in zero field and the dashed line is for the t = 0.2 
junction when in the presence of a magnetic field that is aligned with the vector order pa- 
rameter A. If one considers a magnetic field perpendicular to the order parameter (h -L A), 
it has no effect on the I-V characteristic. On the other hand, in Fig. HI we display typical 
curves for junctions in which both sides are conventional spin-singlet superconductors and 
junctions connecting a spin-singlet and spin-triplet superconductor. These curves address 
the situation of a potential STM experiment carried out using a superconducting tip made 
out of a conventional superconductor used to probe a superconducting phase of unknown 
symmetry. We remind the reader that, in the case of conventional superconductors (dotted 
lines) and when orbital effects can be ignored, the I-V is not sensitive to applied magnetic 
fields. On the other hand, the solid lines correspond to singlet-triplet junctions, that are 
insensitive to the orientation of the vector order parameter on the triple-pairing side of the 
junctions, and their current amplitude is found to be systematically smaller than in the case 
of the respective singlet-singlet junctions. The 'sub-gap' structure shows only two steps and 
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Fig. 4. I-V characteristics of S-T junctions with and without magnetic field (dashed and solid lines 
respectively, h = 0.2). The lower solid and dashed curves are for t = 0.2 (without and with 
magnetic field) and the upper solid curve is for t = 0.5. The dotted lines are (i) the straight 
unitary slope line given for reference and (ii) the S-S characteristics for similar-parameters 
junctions (upper t = 0.5 and lower t = 0.2) given for comparison purposes. 

the current is zero when eV < /^Triplet (the vector order parameter on the spin-triplet side 
of the junction). Concerning the effects of an applied magnetic field, the curves remain 
unchanged if the field is applied parallel to the direction of the vector-order-parameter, but 
show instead a Zeeman effect if the field is perpendicular to it (dashed line). 

Regarding the organic superconductors, 1 - 1 the experiments show that for magnetic fields 
along the direction of the conducting chains (a) the upper critical field is possibly paramag- 
netically limited for small fields (before crossing the upper critical field along b' that is never 
paramagnetically limited). In that range we could assume that the direction of the order 
parameter is fixed respect to the lattice and does not follow the applied field. 38 -* With this 
geometry, a Zeeman splitting of the differential conductance peak should be observed in a 
normal-tip STM experiment. As the field is rotated the splitting would be suppressed and for 
a magnetic field oriented parallel to the b' crystalline-axis there should be no Zeeman effect. 
The disappearance of splitting, even as the field is possibly being increased, would constitute 
a clear signature of spin-triplet superconductivity. Similarly, an s-wave-tip STM would also 
be a direct probe for spin-triplet order. When a field is applied along the b' crystalline-axis, 
a Zeeman splitting would occur. This would constitute a clear sign of unconventional su- 
perconductivity since such an effect does not take place for standard BCS superconductors. 
The b' direction is the one on which the upper critical field is not paramagnetically limited, 
so relatively large fields could be applied in order to obtain a clear signal, and as the field 
alignment changes the splitting should disappear. No successful attempts of this kind of 
experiments were as yet made in the case of the quasi-one- dimensional organic salts, but 
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efforts in this direction are on their way (recently, preliminary experiments involving junc- 
tions between two Bechgaard salts were performed, and they showed a number of puzzling 
features including a zero-bias conductance peak and zero excess current). 55 - 1 We certainly 
hope that further tunnelling experiments along the lines described on these notes can be 
performed. 

§5. Conclusions and perspectives 

In these notes we have presented a new formalism, based on Keldysh technique, to 
compute the full current-voltage characteristics of point contact junctions between normal 
or superconducting (singlet or triplet) leads. This formalism is flexible enough to easily 
incorporate the effects of magnetic fields or temperature. As a possible application of this 
formalism we have shown how it can be used to devise a tunnelling experiment allowing to 
probe for the triplet nature of the superconducting phase in the Bechgaard salts. 

Clearly, this formalism can be extended in many ways. One of the most interesting 
extensions, in view of recent tunnelling experiments, 55 -* is to take into account a finite region 
for the tunnelling contact, an intermediate between the point contact and the planar junction. 
Other extension, needed for example in the case of ruthenates, would be to incorporate the 
modifications of the superconducting order parameter close to the surface. Such a difficulty 
could be avoided for the organics given the large mass anisotropy, but is more crucial for an 
isotropic two dimensional system, since surfaces are pair breaking for triplet superconductors 
when the momentum is perpendicular to the surface. Finally, extensions to other pairing 
symmetries such as rf-wave are also possible; such extensions are, of course, exciting questions 
to be covered in future studies. 
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